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Abstract 

The propagation of massive coloured particles in stochastic background chromo- 
electric field is studied using the semiclassical equations of motion. Depending on the 
nature of the stochastic background we obtain the formulae for the energy losses of 
heavy coloured projectile in nonperturbative hadronic medium and for the energy gains 
in the stochastic field present, e.g., in the turbulent plasma. The result appears to be 
significantly dependent on the form of the correlation function of stochastic external 
field. 
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1. It has become a part of the standard QCD wisdom that the physics of heavy quarks 
provides unique possibilities for studying both perturbative and non-perturbative phenom- 
ena. From the point of view of perturbation theory, due to the asymptotic freedom the 
processes involving heavy quarks tend to evolve at small distances [1]. At the same time the 
heavy quarkonia are a unique testing device for studying the nonperturbative effects. One of 
the best examples is the appearence of the vacuum condensates in the description of heavy 
quarkonia physics in the framework of sum rules method [2]. 

Of a particular interest to high energy physics are the unique abilities of heavy quarkonia 
to probe the properties of a complicated medium arising in the high energy nucleus-nucleus 
and hadron-nucleus collisions. In particular charmonium, when propagating through the 
quark-gluon plasma, effectively measures its temperature giving a possibility of its determi- 
nation through the J/^ to J/\I/' ratio [3]. Another sensitive point, which was a subject of 
recent discussions, is a drastic difference between the interaction of heavy quarkonium with 
the medium, when the heavy quarkonium is in the singlet and octet colour states respectively 
(see e.g. [4]). In particular, such effects can be very important for the description of the hard 
region of spectra in the longitudinal momenta of charmonium, produced in hadron-nucleus 
collisions [5]. Roughly speaking, the small dipole interaction of the singlet projectile with the 
coloured medium is supposed to be responsible for colour transparency, and that of an octet 
one for colour opacity. Let us note, that the production of an initial octet then emitting 
a gluon and becoming a singlet constitutes a part of the perturbative contribution to the 
production of the heavy vector mesons [6]. 

Another situation of considerable interest to the high energy physics is a propagation of 
a heavy quark through the hadronic (nuclear) medium. As the process of creation of heavy 
quarks can often be reliably calculated in perturbative QCD, one hopes to get an information 
on the medium through which the heavy quark (heavy quark jet) propagates. 

One of the models for the interior of the hadronic medium is that of the soft stochastic 
colour fields giving rise to the phenomenologically introduced condensates. Different variants 
of this model were successfully applied both for the explanation of the origin of confinement 
[7] and for more practical calculations such as those of the meson spectra (see e.g. [8,9]) and of 
the soft hadron-hadron scattering cross-sections [10]. Below we shall consider the application 
of this language to a case of a massive colour projectile in a stochastic chromoelectric field. 
In this situation it is interesting to know, how the interaction with the medium changes 
the properties of the incident particle (for example, its energy, etc.). We shall consider the 
changes of the projectile energy due to the interaction with the stochastic chromoelectric 
external field. It will be shown that the particle can either lose or gain energy depending 
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on the model adopted for the description of the hadronic medium. The mechanism of 
this stochastic decceleration (acceleration) is completely analogous to the phenomenon of 
stochastic plasma heating well-known in plasma physics (see, e.g., [11]). In this case when the 
stochastic electric fields are present in the plasma (for example, this is a case for the turbulent 
plasma), the charged particles are accelerated by it. This model was succsefully applied to 
a description of the cosmic rays spectrum, where the charged particles are accelerated by a 
turbulent cosmic plasma (see, e.g, [12]). 

2. Let us consider the propagation of a particle having a non-abelian charge in some 
representation of a gauge group in the external non-abelian field. The cooresponding quasi- 
classical equations of motion are well-known [13,14] and read 

mil, = gQ a (r)G%(x(r))x^ 
Q a (r) = -gf abc x^r)A^x(r))Q c (r) 

where the derivative in the left hand side is taken with respect to a proper time r, m is a 
mass of an incident particle, g is a colour charge, Q a are the components of a colour spin, f a b c 
are the structure constants of a corresponding Lie algebra, x M (r) is the particle trajectory 
in the space-time, w M (r) is a corresponding four- velocity, A h ^ are the external gauge field 
potentials forming a stochastic ensemble and G a is a corresponding field strength. The first 
equation describes the change of energy-momentum of the particle due to the interaction 
with the external field, and the second one describes the precession of the colour spin in the 
external nonabelian field. 

In the following we shall confine our consideration to the one-dimensional case and look 
at the changes in the energy of the massive coloured projectile. 

It is convenient to rewrite the system of equations (1) in the laboratory system: 

= gCr(t)EZ(t,x)(l-vl)W 

^ = gQ a (t)E a x (t,x)v x 

^ = gf abc (v x A b x (t,x)-A b (t,x))Q c (t), (2) 

where the particle is moving along the x-axis with an instantoneous velocity v x (t), E x (t, x) is 
a chromoelectric field strength, £ is energy of the particle. The calculations presented below 
are straightforwardly generalizing those by Sturrock [11] in his discussion of a stochastic 
plasma heating on a relativistic non-abelian case. 

3. Let us now discuss the description of the ensemble of the stochastic chromoelectric 
field. We assume it is a stationary ensemble characterized by the lowest order correlators for 
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the group SU(N C ): 



(E a x (t,x)) 
(E a x (t,x)E b x (t',x>)) 







1 1 
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)S ab R{t-t',x-x'), 



(3) 



where the brackets denote the averaging over the stochastic ensemble. The dimension- 
ful quantity (E 2 ) is an average value of the chromoelectric field squared (nonperturba- 
tive condensate). For the dimensionless correlation function R one has to introduce some 
parametrization. The Eucledian counterpart of this correlation was measured on the lattice 
[15]. Different parametrizations of it were used in phenomenological applications (I.V.Andreev 



Of a crucial importance for the following discussion is the sign of the condensate (E 2 ). 
In the (by now) standard approach to the description of the nonperturbative phenomena in 
terms of condensates this quantity is negative. This shows itself, for example, in the unusual 
sign of the quadratic Stark effect due to the stochastic chromoelectric vacuum fields in heavy 
quarkonia [16]. At the same time the sources of chromoelectric fields like a system of heavy 
charges considered in the discussions of the QCD analog of the Landau-Pomeranchuk effect 
[17] or the instabilities of the quark-gluon plasma correspond to a positive sign of (E 2 ). Below 
it will be shown, that it is precisely this difference in sign that leads to stochastic decceleration 
or acceleration of the projectile particle in the external stochastic chromoelectric field. 

We conclude this section by some technical remarks. In the following we shall work in the 
lowest (second) order in the chromoelectric field and suppose that the autocorrelation time 
of the stochastic external field is less than a characteristic time scale of a problem (i.e. r c < 1 
fm). Another point is that as the third equation in (3) contains the gauge potentials it is 
necessary to express them through the electric field strength. The most straightforward way 
to do it is to use a Fock-Schwinger gauge x^A" = 0, in which the potentials are expressed 
through the field strength by the formula 



which we shall use below. 

4. Let us now look at the evolution of the average particle energy {£). From the Eq. (2) 
we see that in order to calculate the right hand side of the second equation in the second 
order in E a (t,x) one has to calculate the changes in the coordinate, charge and velocity in 



[6], [8]). 
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Then the equation for the charge evolution takes the form 
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the first order in the chromoelectric field. Assuming the unperturbed motion of the form 

x o(t) — v ot, one gets 

x(t) = v t + ^-Q a Q (l - v f 2 f dt'{t - t')E a x (t', vol/) 
m Jo 

v x (t) = vo + ^-(l-vo) 3/2 J*dt'E a x (v t',t') 

Q a (t) = Q a + ^f abc I' daa f dt't'E b x (t',v t')Q c (6) 
m Jo Jo 

where Qq is an initial colour spin. Substituting the formulas (6) into Eq.(2) and averaging 
over the stochastic external field, we get in the leading order in E a : 

<§> = s^ry «' - *"W f - <'>• *<* - '» 

+(1 - ^ 2 ) 3/2 (Qo) 2 / dt'v (t - t')R'((t - t'),v (t - 0) 
Jo 

rl roo 

+ 2v f ab cS ab Q c / daa / dt't'R(t - t', v (t - t'))} (7) 

JO JO 

where the derivative in the second term is taken with respect to a spatial coordinate and 
the infinite limit for the integration over time follows straightforwardly from the assumption 
that the autocorrelation time of the external field is small on the characteristic time scale of 
a problem. 

From the antisymmetry of the structure constants f a b c it immediately follows that the 
third contribution is equal to zero and we finally get 

(f > = 3^ E2) 1) (1 - ^o 2 ) 3/2 (Qo) 2 f dt(R(t, v t) + v tR'(t, v t)) (8) 

The formula (8) clearly shows that the sign of the energy evolution is determined by that 
of an average square of the stochastic chromoelectric field (condensate) (E 2 ). For the usual 
nonperturbative condensate (negative (E 2 )) we have energy losses (stochastic decceleration), 
and in the opposite case of a positive (E 2 ) we have the standard stochastic acceleration well- 
known in plasma physics. 

In order to study the dependence of the energy loss (gain) of the colour projectile on the 
autocorrelation time of the stochastic chromoelectric field t c we have to introduce a specific 
parametrization of the correlation function R(t—t', x—x'). Below we shall see that a choice of 
this parametrization can dramatically change the resulting expression for the energy losses. 

Let us consider two possible parametrizations of the correlation function R: 

R l = exp{-{{t - t'f + (x - x'f )/tl) (9) 
R 2 = eX p(-(( t -t') 2 -(x-x') 2 )/t 2 c ) (10) 
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Let us stress that the second Lorentz-invariant expression provides a natural analytical con- 
tinuation of one of the parametrizations for the correlation function studied in Eucledian 
formalism. In Eq. (9) we have the correlation function R taken on the nonperturbed particle 
trajectory. Thus we get 

R^t) = exp(-(l + v 2 )t 2 /t 2 c ) (11) 
R2(t) = exp(-(l-v 2 )t 2 /t 2 c ) (12) 

Substituting these parametrization into Eq. (9) we get 

d£ 2-n 3 / 2 a s {E 2 )r l-v 2 /2 Q 2 
{ dt h ~ 3m [ l+v 2 ' N*-l 

and 

dS _ 2^ 2 a s {E 2 )r c Q 2 

{ W 2 ~ 3^ Nf^l (14) 

We can conclude that the usual description of the hadronic medium (negative (E 2 )) leads to 
the energy loss of a colour projectile in the stochastic colour medium thus softening the energy 
spectrum of the produced hadrons. At the same time the stochastic acceleration of heavy 
quarkonium by stochastic fields in quark-gluion plasma originating from its instabilities or 
even turbulence will harden this spectrum. 

Denoting the initial energy of the projectile by E and recalling, that the conserved square 
of the colour spin is equal to [14] 

(Q a ) 2 = 3C A 

where Ca is a Casimir operator of the corresponding representation of the colour qroup 
(for SU(N C ) Cy = N c , Cf = (N 2 — 1)/2N C for the vector and fundamental representations 
respectively) we obtain the final formulas 

d£ 2^ 2 a s (E 2 }r C r 2E 2 /2 C A 

m Nf^l (15) 

and 

d£. 2n 3 / 2 a s (E 2 )r c C A 

{ W 2 ~ m Nf^l (16) 

We see that the energy dependence of the energy losses rate has a dramatic dependence on 
the form of the correlation function. 

For the intuitively appealing (although Lorentz-noninvariant) choice Ri the rate of the 
energy losses (gains) decays (at sufficiently high energies of the incident particle Eq) as Eq 3 
and thus is numerically negligible. 
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The second Lorentz-invariant choice R2 leads to an energy independent energy loss (gain) 
rate and thus can be quantitatively important. 

Let us estimate the corresponding energy losses for the charmed quark (fundamental 
representation of SU(3)) and the J/^f meson in the octet state (adjoint representation of 
577(3)) in the case of a constant energy change rate. Using 

(a s E 2 ) = -\(a s G 2 ), 

m c = 1.5 GeV , (a s G 2 ) = 0.042 GeV 4 and t c = 0.3 fm we get 

A =98^ 
{ dt )q ™ fm 

and 

,d£ s MeV 

= 110 ^r 

Taking into account the additional numreical uncertainities (the value of t c , etc) we see, 
that the nonperturbative decceleration of the described type can be around 10-20 percent 
with respect to a naively expected stringy energy loss in the hadronic medium. 

Concluding this section we note, that the formulas (15) and (16) illustrate (within log- 
arithmic accuracy) the borderline energy dependences following from the simplest choice of 
a correlation function characterizing the ensemble of the stochastic chromoelectric external 
fields. 

5. Let us briefly summarize the results. Considering the quasiclassical equations of 
motion for a heavy coloured particle propagating in the external stochastic chromoelectric 
field we have calculated the energy evolution rate. The calculations are generalizing to a 
non-abelian relativistic case those explaining the stochastic heating phenomenon in plasma 
physics. We have shown that the cases of a projectile decceleration and acceleration are dis- 
tinguished by the sign of an chromoelectric field condensate. The negative sign corresponds 
to stochastic decceleration (energy loss) and the positive one to stochastic acceleration (en- 
ergy gain). The energy dependence of the rate of the energy loss turns out to be very sencitive 
to the form of the correlation function of the stochastic external chromoelectric field. 

The first possibility (negative (E 2 )) is believed to be realized in the nonperturbative 
colour medium leading for example to an unusual sign of a level shift of heavy quarkonium 
due to a quadratic Stark effect [16]. The experimental prediction will be the corresponding 
softening of the energy spectrum of the produced hadrons. 

The most natural situation when the second possibility (positive (E 2 )) is realized is for 
hadron production in the quark-gluon plasma in the presence of an instability (turbulence) 
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or a propagation of coloured projectile originating from the hard stage of a collision through 
it. The experimental prediction is the hardening of the energy spectrum of the produced 
hadrons. This phenomenon is analogeous to the stochastic acceleration of charged particles 
in the cosmic plasma responcible for the formation of the energy spectrum of cosmic rays 
[12]. 
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